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We report the measurement of v-e elastic scattering from 8 B solar neutrinos with 3 MeV energy 
threshold by the Borexino detector in Gran Sasso (Italy). The rate of solar neutrino-induced electron 
scattering events above this energy in Borexino is 0.217 ± 0.038(stai) ± 0.008(syst) cpd/lOOt, which 
corresponds to "i&gj = 2.4 ± 0.4± O.lxlO 6 cm" 2 s -1 , in good agreement with measurements from 
SNO and SuperKamiokaNDE. Assuming the 8 B neutrino flux predicted by the high metallicity Standard 
Solar Model, the average 8 B v e survival probability above 3 MeV is measured to be 0.29±0.10. The 
survival probabilities for 7 Be and 8 B neutrinos as measured by Borexino differ by 1.9 a. These results are 
consistent with the prediction of the MSW-LMA solution of a transition in the solar v e survival probability 
Pee between the low energy vacuum-driven and the high-energy matter-enhanced solar neutrino oscillation 
regimes. 



PACS numbers: 14.50. St, 26.65.+t, 95.55.Vj, 29.40.Mc 

INTRODUCTION 

Solar 8 B-neutrino spectroscopy has been so far per- 
formed by the water Cerenkov detectors KamiokaNDE, Su- 
perKamiokaNDE, and SNO Q-Q. The first two experi- 
ments used elastic v-e scattering for the detection of neu- 
trinos, whereas SNO also exploited nuclear reaction chan- 
nels on deuterium with heavy water as target. These experi- 
ments provided robust spectral measurements with ~5 MeV 
threshold or higher for scattered electrons; a recent SNO 
analysis reached a 3.5 MeV threshold [5]. 

We report the first observation of solar 8 B-neutrinos with 



a liquid scintillator detector, performed by the Borexino ex- 
periment 0, 0] via elastic v-e scattering. Borexino is the 
first experiment to succeed in suppressing all major back- 
grounds, above the 2.614 MeV 7 from the decay of 208 TI, 
to a rate below that of electron scatterings from solar neu- 
trinos. This allows to reduce the energy threshold for scat- 
tered electrons by 8 B solar neutrinos to 3 MeV, the lowest 
ever reported for the electron scattering channel. To facil- 
itate a comparison to the results of SuperKamiokaNDE 0] 
and SNO D2O phase [4], we also report the measured 8 B 
neutrino interaction rate with 5 MeV threshold. 

Since Borexino also detected low energy solar 'Be neutri- 
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nos [8J, |9), this is the first experiment where both branches 
of the solar pp-cycle have been measured simultaneously in 
the same target. The large mixing angle solution (LMA) 
of the MSW effect fiol ] predicts a transition in the v c sur- 
vival probability from the vacuum oscillation regime at low 
energies to the matter dominated regime at high energies. 
Results on solar 7 Be and 8 B neutrinos from Borexino, com- 
bined with prediction on the absolute neutrino fluxes from 
the Standard Solar Model confirm that our data 

are in agreement with the MSW-LMA prediction within la. 



EXPERIMENTAL APPARATUS AND ENERGY 
THRESHOLD 

The Borexino detector is located at the underground Lab- 
oratori Nazionali del Gran Sasso (LNGS) in central Italy, at a 
depth of 3600 m.w.e.. Solar neutrinos are detected in Borex- 
ino exclusively via elastic v-e scattering in a liquid scintil- 
lator. The active target consists of 278 1 of pseudocumene 
(PC, 1,2,4-trimethylbenzene), doped with 1.5 g/l of PPO 
(2,5-diphenyloxazole, a fluorescent dye). The scintillator is 
contained in a thin (125 jim) nylon vessel of 4.25 m nominal 
radius, and is shielded by two concentric PC buffers (323 
and 567 1) doped with 5.0 g/l of a scintillation light quencher 
(dimethylphthalate). Scintillator and buffers are contained 
in a Stainless Steel Sphere (SSS) with a diameter of 13.7 m 
and the scintillation light is detected via 2212 8" photo- 
multiplier tubes (PMTs) uniformly distributed on the inner 
surface of the SSS. The two PC buffers are separated by 
a second thin nylon membrane to prevent diffusion of the 
radon emanated by the PMTs and by the stainless steel of 
the sphere into the scintillator. The SSS is enclosed in a 
18.0 m diameter, 16.9 m high domed Water Tank (WT), 
containing 2100t of ultra-pure water as an additional shield- 
ing against external gamma- and neutron background. 208 
8" PMTs in the WT detect the Cerenkov light produced by 
muons in the water shield, serving as a highly efficient muon 
veto. A complete description of the Borexino detector can 
be found in Ref. [g]. 

Scintillator detectors, with their high light yield, are sen- 
sitive to lower energy events than Cerenkov detectors. In 
this analysis, the 3 MeV energy threshold is imposed mainly 
by the 2.614 MeV 7 -rays from the decay of 208 TI ( 232 Th 
chain, Q=5.001MeV) in the PMTs and in the SSS and 
by the finite energy resolution of the detector: a tail of 
2.6 MeV 7 events leaks at higher energies, along with a 
very small percentage of combined 7's from 214 Bi ( 238 U 
chain, Q=3.272 MeV). The 3 MeV energy threshold elimi- 
nates sample contamination from such events. 

Potential background sources above 3 MeV include the 
radioactive decays of residual 214 Bi and 208 TI within the liq- 
uid scintillator, decays of cosmogenic isotopes (see Table [I] 
later), high energy 7-rays from neutron capture, and cosmic 
muons. No background from a decays is expected at these 
energies, since the light quenching of as in organic liq- 



uid scintillators reduces their visible energy in the electron- 
equivalent scale below 1 MeV. A measurement of 8 B neutri- 
nos with a 3 MeV energy threshold is contingent upon high 
radiopurity of the scintillator target. 238 U and 232 Th con- 
centrations in the Borexino scintillator have been measured 
at (1.6±0.1)xl0- 17 g/gand (6.8±1.5)x 10~ 18 g/g, respec- 
tively, and record low levels of backgrounds in the energy 
range 0.2-5.0 MeV have been reported [fil, E|- The dom- 
inant background in the energy range of interest for solar 
8 B neutrinos originates from spallation processes of high en- 
ergy cosmic muons. This paper demonstrates that, thanks 
to the LNGS depth and the Borexino muon veto system, 
cosmogenic background can be reduced below the rate of 
interaction of 8 B neutrinos, thus allowing the neutrino rate 
to be measured. 



ENERGY AND POSITION RESPONSE OF THE 
DETECTOR AND ASSOCIATED SYSTEMATICS 

We tuned the response of energy and position reconstruc- 
tion algorithms with a dedicated calibration campaign. We 
used an off-axis source insertion system designed to posi- 
tion radioactive and/or luminous sources at several locations 
throughout the detector active target. The source position 
can be measured with a set of stereoscopic cameras installed 
on the SSS with an uncertainty of ±2 cm in x, y, and z 14 1 
(±3. 5cm in radius). 



As explained above, Borexino has a unique sensitivity to 
electron scattering from low energy solar neutrinos, thanks 
to its unmatched record on background below the natural 
radioactivity barrier. To preserve this capability, the off-axis 
source insertion system was designed to respect stringent 
limits on leak tightness and cleanliness of the mechanics in 
contact with the liquid scintillator. A detailed description of 
the technique used for the calibration and for the position 
reconstruction in Borexino is in preparation. 



Energy Scale 

Calibration of the energy scale allows to establish with 
high confidence the energy threshold for the 8 B neutrino 
analysis and the error in its determination, and to calibrate 
the energy scale to allow the energy spectrum of the elec- 
trons scattered by 8 B neutrinos to be determined. Distor- 
tions of the energy scale are due to physical effects (quench- 
ing), geometrical effects (light collection), and to the elec- 
tronics, which was designed for optimal performance in the 
low-energy range of 7 Be neutrinos, in a regime where a 
single photoelectron is expected for each PMT. At higher 
energies, the electronics response to multiple photoelectron 
hits on a single channel is not linear. For each triggered 
channel, the charge from photoelectrons in a 80 ns gate is 
integrated and recorded, but photoelectrons in the follow- 
ing 65 ns dead time window are lost. The resulting fraction 
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of lost charge increases with energy and can reach ~10% 
at ^10 MeV. Moreover, the number of detected photoelec- 
trons depends on the event position in the active volume, 
due to differences in PMT coverage. For an accurate deter- 
mination of the energy scale, the dominant non-linearities 
have been reproduced with a Monte Carlo simulation and 
with calibration measurements. 

To avoid contamination during calibration, j3 sources were 
not put in direct contact with the scintillator. Instead, we 
calibrated the response of the detector with encapsulated 
7 sources. The scintillation induced by 7-rays is due to 
the ionizing tracks of the secondary electrons, thus, the two 
energy scales are closely related. 

Establishing a correlation between the (3 and 7 energy 
scales required extensive simulations with the G4Bx Monte 
Carlo code. G4Bx is based on Geant4 [lil, 16| and simu- 
lates in detail all of the detector component, and includes 
scintillation, Cerenkov photon production, absorption and 
scattering of light in the scintillator and in the buffer, as 
well as the PMT response. Each secondary electron in a 
7-induced Compton electron cascade is affected by energy- 
dependent ionization quenching, which amplifies the distor- 
tion in the 7 energy scale. The quenching effect is modeled 
with the Birks formalism [17]. A second package, BxElec, 
simulates in detail the response of the electronics. Finally, 
Monte Carlo data are processed by the same reconstruction 
code used for real data. A detailed description of the Monte 
Carlo codes and of the energy reconstruction algorithm is in 
preparation. 

To calibrate the detector energy response to 8 B neutri- 
nos, we used an 241 Am 9 Be neutron source positioned at the 
center of the detector and at several positions at 3 m radius. 
Neutron capture on 1 H and on 12 C in the scintillator results 
in the emission of 7-rays from the 2.223 MeV and 4.945 MeV 
excited states, respectively. In addition, neutron capture on 
the stainless steel of the insertion system produces 7-rays 
from the 7.631 MeV ( 56 Fe) and 9.298 MeV ( 54 Fe) excited 
states. We validate the Monte Carlo code by simulating the 
four 7-rays in both the positions. In Figure [1] we show the 
results of the calibration of the 7-equivalent energy scale 
in the detector center. Monte Carlo simulations reproduce 
7 peak positions and resolutions at o\ — 1% precision in 
the detector center (as shown in Figure [1]), and at 02 — 
4% precision at 3 m from the detector's center. Assuming 
the same accuracy for the /3-equivalent energy scale, we ex- 
trapolate it by simulating electrons uniformly distributed in 
the scintillator, and then selecting those with reconstructed 
position within the fiducial volume. The error on the en- 
ergy scale is obtained with a linear interpolation from o\ in 
the detector center to 02 at 3 m, along the radius. The j3- 
equivalent energy scale, in the energy region above 2 MeV, 
can be parametrized as: 

N = a-E + b, (1) 
where N is the number of photoelectrons (p.e.) detected by 
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Figure 1: Black dots are the measured peak positions of 7 
radiation induced by neutron captures in 1 H (2.223 MeV), 
12 C (4.945 MeV), 56 Fe (7.631 MeV) and 54 Fe (9.298 MeV) in 
the detector center. Red line is the Monte Carlo prediction 
for 7 rays generated in the detector center. 



the PMTs, a=459±llp.e./MeV and 6=115±38p.e.. The 
non-zero intercept b is related to the fact that this descrip- 
tion is valid only in this energy range and that the overall 
relation between N and E is non linear. 

The anticipated 3 MeV (5 MeV) energy threshold for the 
8 B analysis corresponds to 1494 p.e. (2413 p.e.) within a 
3 m radial distance from the center of the detector. The un- 
certainty associated to the 3 MeV (5 MeV) energy threshold 
is obtained by propagating the errors of Eq.[TJand is equal 
to 51 p.e. (68 p.e.). 



Vertex Reconstruction 

The positions of scintillation events are reconstructed 
with a photon time-of-flight method. We computed with 
G4Bx a probability density function (PDF) for the time of 
transit of photons from their emission point to their detec- 
tion as photoelectron signals in the electronics chain. We 
refined the PDF with data collected in the calibration cam- 
paign. Event coordinates (xq, yo, zq) and time (to) are 
obtained by minimizing: 



C(x 0) yo, z , t ) = Y[ PDF 



to 



do.i ■ n c ff 



(2) 



where the index i runs over the triggered PMTs, ti is the 
time of arrival of the photoelectron on the i th electronic 
channel, and drj,j is the distance from the event position 
and the i th PMT. n c g is an empirically-determined effec- 
tive index of refraction to account for any other effect that 
is not accounted for in the reconstruction algorithm but im- 
pacts the distribution of PMT hit times, both in the optics 
(e.g. Rayleigh scattering) and the electronics (e.g. multiple 
photoelectron occupancy). 
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Figure 2: Ratio of the reconstructed radial position of 7 events 
from the 241 Am 9 Be source in Borexino to the source radial 
position measured by the CCD camera system, as a function 
of the measured charge. 



8 B-NEUTRINO FLUX 

We report our results for the rate of electron scattering 
above 3MeV from 8 B neutrino interactions in the active 
target. We also report the result above the threshold of 
5 MeV, to facilitate the comparison with results reported by 
SNO 0| and SuperKamiokaNDE phase-l @] at the same 
threshold. This energy range is unaffected by the scintil- 
lator intrinsic background, since the light quenching effect 
reduces the visible energy of 208 TI (Q=5.001 MeV) from 
232 Th contamination in the scintillator below the energy 
threshold of 5 MeV. 

The analysis in this paper is based on 488 live days of data 
acquisition, between July 15, 2007 and August 23, 2009, 
with a target mass of 100 1, defined by a fiducial volume 
cut of radius 3 m. The total exposure, after applying all the 
analysis cuts listed in the next section, is 345.3 days. 



The Borexino electronics records the time of each de- 
tected photoelectron introducing a dead time of 145 ns 
after each hit for each individual channel. Therefore, the 
timing distribution is biased at high energy, where multiple 
photoelectrons are detected by each channel, and the po- 
sition reconstruction is energy dependent. To measure this 
effect, we deployed the 241 Am 9 Be neutron source at the six 
cardinal points of the sphere defining the fiducial volume, 
i.e. those points lying on axis through the center of the 
detector, with off-center coordinates from the set i=±3m, 
y=±3m, and z=±3m. The recoiled proton from neutron 
scattering allows us to study the reconstructed position as 
function of the collected charge up to ~5000 p.e.. Figure[2] 
shows the ratio of measured versus nominal position of the 
241 Am 9 Be source. This data was used to define the fiducial 
volume i? nom <3m. The non-homogeneous distribution of 
live PMTs, in particular the large deficit of live PMTs in 
the bottom hemisphere is responsible for the different 
spatial response at mirrored positions about the x-y plane. 
Thus, as shown in Figure [2] two radial functions have been 
defined for positive and negative z positions. 



After all post-calibration improvements to the event re- 
construction algorithm, typical resolution in the event po- 
sition reconstruction is 13±2cm in x and y, and 14±2cm 
in z at the relatively high 214 Bi energies. The spatial reso- 
lution is expected to scale as 1/VN where N is the num- 
ber of triggered PMTs, and this was confirmed by deter- 
mining the 14 C spatial resolution to be 41±6cm (la) at 
140 keV (HQ. Systematic deviations of reconstructed posi- 
tions from the nominal source position are due to the 3.5 cm 
accuracy of the CCD cameras in the determination of the 
calibration source position, and in 1.6 cm introduced by 
the energy dependency. The overall systematics are within 
3.8cm throughout the 3m-radius fiducial volume. 



Muon Rejection 

The cosmic muon rate at LNGS is 1.16±0.03 m _2 hr _1 
with an average energy of 320±4 sta t± Hsys GeV [18]. Each 
day, ^4300 muons deposit energy in Borexino's inner detec- 
tor. Depending on deposited energy and track length, there 
is a small but non-zero chance that a cosmic muon induces 
a number of photoelectrons comparable to the multi-MeV 
electron scatterings of interest for this analysis, and is mis- 
taken for a point-like scintillation event. A measurement of 
the neutrino interaction rate in Borexino requires high per- 
formance rejection of muon events and an accurate estimate 
of the muon tagging efficiency. 

As mentioned earlier, the Borexino WT is instrumented 
with 208 PMTs to serve as a muon veto. If an Inner De- 
tector (ID) event coincides in time with an Outer Detector 
(OD) trigger (i.e. more than 6 PMTs in the WT are hit 
within a 150 ns window), the event is tagged as muon and 
rejected. However, the OD efficiency is not unity and de- 
pends on the direction of the incoming cosmic muon. 

In addition, we perform pulse-shape discrimination on the 
hit time distribution of inner detector PMTs, since for track- 
like events, like muons, such distribution generally extends 
to longer times than for point-like events, like /3-decays 
and v — e scattering. We exclude muons from the event 
sample in the energy range of interest (3.0-16.3 MeV, or 
1413-6743 p.e.) by imposing the following requirements (ID 
cuts): 

• The peak of the reconstructed hits time distribution, 
with respect to the first hit, is between ns and 30 ns. 

• The mean value of the reconstructed hits time distri- 
bution, with respect to the first hit, is between ns 
and 100 ns. 

The efficiency of the selection cuts was evaluated on a 
sample of 2,170,207 events, identified by the OD as muons. 
Only 22 of these events, a fraction of (1.0±0.2)xl0~ 5 , sur- 
vive the ID cuts in the energy and spatial region of interest, 
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and are tagged as possible scintillation events. 

We do not have an absolute value for the OD muon 
veto efficiency, but we estimate it to be larger than 
99%, from G4Bx simulations. The residual muon rate, 
due to the combined inefficiency of the two tagging sys- 
tems, taking into account the fact that the two detectors 
are independent, is (4.5±0.9)x 10~ 4 muons/day/lOOt, or 
(3.5±0.8)xl0~ 4 muons/day/lOOt above 5 MeV. 

Cosmogenic Background Rejection 

Fast cosmogenic veto 

Table [T] presents a list of expected cosmogenic isotopes 
produced by muons in Borexino. The short-lived cosmo- 
genics (r < 2s), as well as the 7-ray capture on 12 C, are 
rejected by a 6.5s cut after each muon, with a 29.2% frac- 
tional dead time. Figure [3] shows the time distribution of 
events following a muon. The data is well fit by three ex- 
ponentials with characteristic times of 0.031±0.002s ( 12 B), 
0.25±0.21s ( 8 He, 9 C, 9 Li), 1.01±0.36s ( 8 B, 6 He , 8 Li), in 
good agreement with the lifetimes of the short-lived isotopes 
(see Table[T]). From the fit we estimate the production rates 
of these cosmogenic isotopes in Borexino. We conclude 
that rejection of events in a 6.5 s window following every 
muon crossing the SSS reduces the residual contamination 
of the short lived isotopes to (1.7±0.2)xlCT 3 cpd/100t 
((1.3±0.2)xlCT 3 cpd/100t above 5 MeV). 

The expected rates (R) quoted in Tableware obtained 
by scaling the production rates (R°) measured by Kam- 
LAND with: 

« = «"(!)> (3) 

where E p and <I> M are the Borexino mean muon energy 
(320±4 sta t ± ll S y S GeV) and flux (1.16±0.03 rrr 2 hr" 1 ), 
as measured by MACRO [Hj], and E° (260±4GeV) and i>° 
(5.37±0.41 m" 2 hr _1 ) are the corresponding KamLAND 
values, a is a scaling parameter to relate cosmogenic pro- 
duction rate at different mean energies of the incoming 
muon flux; it is obtained in Ref. [19] by fitting the pro- 
duction yield of each isotope, simulated by FLUKA, as a 
function of muon beam energy. Overall, Borexino data re- 
sults are in agreement with the values quoted in Table [T] 
within 15%. 



Neutron rejection 

The cosmogenic background in Borexino includes decays 
of radioactive isotopes due to spallation processes on the 
12 C nuclei in the scintillator, as well as the 7-rays from the 
capture of neutrons that are common by-products of such 
processes. The capture time for neutrons in the Borexino 




Time [s] 



Figure 3: Cumulative distribution of events with energy > 
3 MeV within a 5 s window after a muon in Borexino. The 
time distribution has been fit to three decay exponentials. 
The ensuing exponential lifetimes are r = 0.031±0.002 s, 
0.25±0.21 s, 1.01±0.36s and corresponds to the contribution 
from 12 B, 8 He + 9 C + 9 Li and 8 B + 6 He + 8 Li, respectively. 
The expected and measured rates for these cosmogenic iso- 
topes are summarized in Table UJ 

scintillator has been measured to be 256.0±0.4/is, using a 
neutron calibration source, and the energy of the dominant 
7-rays from neutron capture on 1 H at 2.223 MeV is below 
the energy threshold of the present analysis. On the other 
hand, the 4.9 MeV 7-rays from neutron captures on 12 C is a 
potential background for this analysis. The rate is estimated 
by scaling the cosmogenic neutron capture rate on 1 H by 
the fraction of captures on 12 C with respect to the total, 
measured with the 241 Am 9 Be neutron source. The neutron 
capture rate on 12 C is 0.86±0.01cpd/100 t. 

The fast cosmogenic veto, described in the Fast Cosmo- 
genic Veto section, rejects neutrons produced in the scin- 
tillator or in the buffer by muon spallation with 99.99% 
efficiency. To reject neutrons produced in water, a second 
2 ms veto is applied after each muon crossing the Water 
Tank only. The rejection efficiency for neutrons produced 
in water is 0.9996. The overall survival neutron rate in 
the energy range of interest and in the fiducial volume is 
(8.6±0.1)xl(r 5 cpd/100 t. 



C identification and subtraction 

A separate treatment is required for long-lived (t>2s) 
cosmogenic isotopes. Since 'Be (t = 76. 9d, Q- 
value=0.9 MeV) and n C (r = 29.4min, Q-value=2.0 MeV) 
are below the energy threshold, we focus on 10 C and 11 Be. 

Taking into account the energy response of Borexino, the 
fraction of the 10 C energy spectrum above 3 MeV is 1.2%. 
When 10 C is produced in association with a neutron, 10 C 
candidates are tagged by the three-fold coincidence with 
the parent muon and subsequent neutron capture in the 
scintillator [20]. The efficiency of the Borexino electronics 
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Isotopes 


T 


Q 


Decay Expected Rate 


Fraction 


Expected Rate > 3 MeV Measured Rate > 3 MeV 






[MeV] 




[cpd/lOOt] 


> 3 MeV 


[cpd/lOOt] 


[cpd/lOOt] 


i '> _ 

B 


0.03 s 


13.4 


P~ 


1.41 ± 0.04 


0.886 


1.25± 0.03 


1.48 ± 0.06 


s He 


0.17s 


10.6 


P~ 


0.026 ± 0.012 


0.898 






C 


0.19s 


16.5 


P + 


0.096 ± 0.031 


0.965 


(1.8 ±0.3 )xl0 _1 


(1.7 ± 0.5)xl0 _1 


9i ■ 

Li 


0.26 s 


13.6 


o— 

P 


n ri7i 4- n fine; 


0.932 








1.11 s 


18.0 


P + 


0.273 ± 0.062 


0.938 






6 He 


1.17s 


3.5 


P~ 


NA 


0.009 


(6.0 ± 0.8)xl0 _1 


(5.1 ± 0.7)xl0 _1 


8 Li 


1.21s 


16.0 


P~ 


0.40 ± 0.07 


0.875 








27.8 s 


3.6 


P + 


0.54 ± 0.04 


0.012 


(6.5±0.5)xl0 _;! 


(6.6±1.8) xl0~ a 


n Be 


19.9 s 


11.5 


P~ 


0.035 ± 0.006 


0.902 


(3.2 ± 0.5)xl0^ 


(3.6±3.5)xl0~- 



Table I: Expected muon-induced contaminants with Q-value > 3 MeV in Borexino. The expected rate is obtained by extrap- 
olating the recently released data by the KamLAND Collaboration [TiJ in accordance with the empirical law defined in Eq 

El 



in detecting at least one neutron soon after a muon has 
been estimated to be 94% by two parallel (1-channel and 
8-channel) DAQ systems that digitize data for 2 ms after 
every OD trigger at 500 MHz. The rate of muons asso- 
ciated with at least 1 neutron, measured by the Borexino 
electronics, is ~67cpd. Thus, to reject 10 C from the anal- 
ysis we exclude all data within a 120s window after a /.t+n 
coincidence and within a 0.8 m distance from the neutron 
capture point. The efficiency of this cut is 0.74±0.11, for a 
0.16% dead time. A time profile analysis of events tagged 
by this veto above 2.0 MeV returns a characteristic time 
of 30±4s, consistent with the lifetime of 10 C, and a to- 
tal 10 C rate of (0.50±0.13)cpd/100t, in production chan- 
nels with neutron emission. Thus, the residual 10 C con- 
tamination from neutron-producing channels above 3 MeV 
is (6.0±0.2)xl0- 3 cpd/100 t. 

The dominant neutron-less 10 C production reaction 
is 12 C(p,i) 10 C. We extrapolated its rate by scaling 
the 12 C(p,d) n C production rate |20|, [2lj|, by the ra- 
tio between the 12 C(p,d) n C and 12 C(p,<) 10 C cross sec- 
tions measured in Ref. [22]. The 12 C(p,i) 10 C rate is 
(0.6±1.8)xl0~ 3 cpd/100 t. The residual background 
above 3 MeV from 10 C is 2.2xl0~ 3 cpd/100 t. 

The overall 10 C rate above 3 MeV, (6.6±1.8) 
xl0 _3 cpd/100 t, agrees with the expected one, quoted 
in TableQ] 

11 Be estimation 

Figure [4] shows the time profile of events within 240 s 
after a muon and within a 2 m distance from its track 
in the entire Borexino active volume (278 t). The effi- 
ciency of the distance cut is assumed to be the same as 
the one measured for cosmogenic 12 B (84%) by perform- 
ing fits to the time distribution of events after a muon 
before and after the track cut. The measured 11 Be rate 
above 3 MeV is (3.6±3.5)x 10~ 2 cpd/lOOt, consistent with 
the (3.2±0.6)xl0~ 2 cpd/100t rate extrapolated from the 
KamLAND measurements (l9| . Since all measured rates de- 
viated less than 18% from the extrapolated value, we adopt 




Time [s] 

Figure 4: Time profile of events with energy > 3 MeV within 
240 s after a muon and within 2 m from its track in the entire 
Borexino active volume (278 t). The distribution has been fit 
to the three decay exponentials as in Figure O plus the 11 Be 
component, with fixed mean-lives. 

the latter as the residual rate of 11 Be in our sample. 



Radioactive Background Rejection 

External background 

The 3 MeV energy threshold is set by the 2.614 MeV 7- 
rays from the /3-decay of 208 TI, due to radioactive contam- 
ination in the PMTs and in the SSS. Above 3 MeV, the 
sources of radioactive background include the radioactive 
decays of residual 214 Bi ( 238 U chain, Q=3.272MeV) and 
208 T! (232 Jh chain> Q =5 .001MeV) in the liquid scintilla- 
tor. The fiducial volume cut is very effective against the 
208 TI and 214 Bi background due to 222 Rn and 220 Rn em- 
anated from the nylon vessel, as well as residual external 
7-ray background. In Figure[5]the radial distribution of all 
scintillation events above 3 MeV has been fit to a model 
which takes into account the three sources of backgrounds: 
a uniform distribution in the detector for internal events, a 
delta-function centered on the vessel radius for the point- 
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Radius [m] 

Figure 5: Fit of the radial distribution for events with 
E>3 MeV. The red line represents the uniformly distributed 
event component in the active mass, the green line the surface 
contamination, and the blue line is external background. 

like radioactive background in the nylon, and an exponential 
for external 7-ray background. All the three components are 
convoluted with the detector response function. From this 
radial analysis we conclude that within the fiducial volume 
there is a small contribution of events from surface con- 
tamination and the exterior of (6.4±0.2)xlCT 3 cpd/100 t 
((3±11) x 1(T 6 cpd/100 1) for events above 3 MeV (5 MeV). 



Bi contamination 

The suppression of 214 Bi from 238 U contamination in 
the FV relies on the 214 F3i- 214 Po delayed coincidence 
(t—237 /j,s). We look for coincidences in time between 20 /js 
and 1.4 ms (e = 0.91) with a spatial separation <1.5 m (e = 
1) and Gatti parameter, a pulse shape discrimination esti- 
mator introduced in Ref. |l,[23j], larger than -0.008 (e = 1) 



Cut 


Counts 


Counts 




3.0-16.3 MeV 5.0 


-16.3 MeV 


All counts 


1932181 


1824858 


Muon and neutron cuts 


6552 


2679 


FV cut 


1329 


970 


Cosmogenic cut 


131 


55 


10 C removal 


128 


55 


214 Bi removal 


119 


55 


208 TI subtraction 


90±13 


55±7 


11 Be subtraction 


79T13 


47±8 


Residual subtraction 


75±13 


46±8 


Final sample 


75±13 


46±8 


BPS09(GS98) 8 B v 


86±10 


43±6 


BPS09(AGS05) 8 B v 


73±7 


36±4 



Table II: Effect of the sequence of cuts on the observed counts. 
The cosmogenic cut introduces a reduction of the detector 
live-time of 29.2%. The resulting effective live-time is 345.3 d. 
The analysis is done with 100 1 fiducial mass target, after the 
FV cut. The expected 8 B counts are calculated from current 
best parameters for the MSW-LMA [11] and Standard Solar 
Models, BPS09(GS98) and BPS09(AGS05) [IM3 



within the fiducial volume throughout the entire data set. 
The 214 Po a-decays are selected in the 0.3-1.2 MeV (e = 1) 
energy range. The remaining contribution of 214 Bi to the v- 
e scattering sample is negligible (l.l±0.4)xl0 _4 cpd/100t. 



Tl contamination 

Amongst the daughters of 232 Th naturally present in the 
scintillator, 208 TI decays are the only ones which contribute 
background above 3 MeV. The parent of 208 TI is 212 F3i a- 
decays into 208 TI with a branching ratio of 36% and a life- 
time of t = 4.47 min. In the second channel with branch- 
ing ratio 64%, 212 Bi /3-decays into 212 Po with a lifetime 
of t = 431 ns. We estimate the 208 TI rate from the fast 
212 Bi- 212 Po coincidences. 212 Bi- 212 Po events are selected 
in a time window between 400 and 1300 ns, with an effi- 
ciency of 0.35, and requiring a maximum spatial distance 
between the two events of 1 m (e— 1). 212 Bi and 212 Po are 
selected in [20-1200] p.e. (e = 1) and [420-580] p.e. (e = 
0.93) energy regions, respectively. The 212 Po-a quenched 
energy is estimated from the 210 Po and 214 Po peaks, opti- 
mal signatures for the a-quenching calibration. We found 
21 212 Bi- 212 Po coincidences in the entire data set, within 
the FV. Accounting for the efficiency of the selection cuts 
and the branching ratios of the 212 Bi decays, this corre- 
sponds to a 208 TI contamination in our neutrino sample of 
29±7 events, or a (8.4±2.0)xl0~ 2 cpd/lOOt rate. 

A summary of the analysis sequence described above is 
shown in TableQT| The energy spectrum of the final sample, 
compared with simulated spectra of 8 B neutrinos and of 
each residual background component listed in Table QTTJ is 
shown in Figure [6] 



NEUTRINO INTERACTION RATES AND ELECTRON 
SCATTERING SPECTRUM 

The mean value for 8 B neutrinos in the sample above 
3 MeV (5 MeV) is 75±13 (46±8) counts. 

The dominant sources of systematic errors are the deter- 
minations of the energy threshold and of the fiducial mass, 



Background 


Rate [10~ 4 cpd/100 t] 






>3 MeV 


>5 MeV 


Muons 


4.5±0.9 


3.5±0.8 


Neutrons 


0.86±0.01 





External background 64±2 


0.03±0.11 


Fast cosmogenic 


17±2 


13±2 


10 C 


22±2 





214 Bj 


1.1±0.4 





208 -p| 


840±20 





"Be 


320±60 


233±44 


Total 


1270±63 


250±44 



Table III: Residual rates of background components after the 
data selection cuts above 3 and 5 MeV. 
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Source 


E>3 MeV 




E>5 MeV 






°"+ 


(7 — 


°+ 


(7 — 


Energy threshold 


3.6% 


3.2% 


6.1% 


4.8% 


Fiducial mass 


3.8% 


3.8% 


3.8% 


3.8% 


Energy resolution 


0.0% 


2.5% 


0.0% 


3.0% 


Total 


5.2% 


5.6% 


7.2% 


6.8% 



Table IV: Systematic errors. 



Data 
MC 

2M TI 
"Be 

External Background 



12 



_l 

14 



Energy [MeV] 

Figure 6: Comparison of the final spectrum after data selec- 
tion (red dots) to Monte Carlo simulations (black line). The 
expected electron recoil spectrum from to oscillated 8 B v in- 
teractions (filled blue histogram), 208 T1 (green), n Be (cyan) 
and external background (violet), are equal to the measured 
values in Tabic HITl 

both already discussed in the previous sections. The first in- 
troduces a systematic uncertainty of +3.6% -3.2% (+6.1% 
-4.8% above 5 MeV). The second systematic source is re- 
sponsible for a ±3.8% uncertainty in the 8 B neutrino rate. 
A secondary source of systematics, related to the effect of 
the energy resolution on the threshold cuts, has been stud- 
ied on a simulated 8 B neutrino spectrum and is responsible 
for a systematic uncertainty of +0.0% -2.5% (+0.0% -3.0% 
above 5 MeV). 

The total systematic errors are shown in Table [TVl 

The resulting count rate with E>3 MeV is: 

0.217 ±0M&(stat)± o Q Q o ° o l(sy st) cpd/lOOt 
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Figure 7: Comparison of the final spectrum after data selec- 
tion and background subtraction (red dots) to Monte Carlo 
simulations (blue) of oscillated 8 B v interactions, with am- 
plitude from the Standard Solar Models BPS09(GS98) (high 
metallicity) and BPS09(AGS05) (low metallicity) , and from 
the MSW-LMA neutrino oscillation model. 

SOLAR 8 B NEUTRINO FLUX AND NEUTRINO 
OSCILLATION PARAMETERS 



The equivalent unoscillated 8 B neutrino flux, derived 
from the electron scattering rate above 5 MeV (Ta- 
ble El), is (2.7±0.4 stat ±0.2 syst )xl0 6 cm-V 1 , in good 
agreement with the SuperKamiokaNDE-l and SNO D2O 
measurements with the same threshold, as reported in 
Table [VT] The corresponding value above 3 MeV, is 
(2.4±0.4 sta t±0.1 sys t)xl0 6 crrr 2 s _1 . The expected value 
for the case of no neutrino oscillations, including the theo- 
retical uncertainty on the 8 B flux from the Standard Solar 
Model [IlHH], is (5.88±0.65)xl0 6 cm^s- 1 and, there- 
fore, solar v e disappearance is confirmed at 4.2o\ 

To define the neutrino electron survival probability P ee 
averaged in the energy range of interest, we define the mea- 
sured recoiled electron rate R, through the convolution: 



dT e f +OC dt(^(t) 

>T C >T J -00 



dt 



between the detector energy response g, assumed gaussian, 
with a resolution depending on the energy, and differential 
rate: 



and with E>5 MeV: 



0.134 ±0m2(stat)ZZ$$(syst) cpd/lOOt 



The final energy spectrum after all cuts and residual back- 
ground is shown in Figure [7] It is in agreement with the 
scenario which combines the high metallicity Standard So- 
lar Model, called BPS09(GS98) [H, and the prediction of 
the MSW-LMA solution. 



dr 

~~pp~\-L ej — 1 

a±e ■' E U >E 



with: 



and: 



dE v (Pee ■ *e + (1 - Pee) ■ ^-r) 



(5) 



* e = ^(E v ,T e ).N e .^-(E v ), (6) 
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3.0-16.3 MeV 


5.0-16.3 MeV 


Rate 

^cxp 
$ ES 


[cpd/100 t] 

[10 6 cm^s- 1 ] 


0.22±0.04±0.01 
2.4±0.4±0.1 
0.88±0.19 


0.13±0.02±0.01 
2.7±0.4±0.2 
1.08±0.23 



Table V: Measured event rates in Borexino and comparison 
with the expected theoretical flux in the BPS09(GS98) MSW- 
LMA scenario [Tot ]. 




Figure 8: Electron neutrino survival probability as function 
of the neutrino energy, evaluated for the 8 B neutrino source 
assuming the BPS09(GS98) Standard Solar Model [IM1 
and the oscillation parameters from the MSW-LMA solution 
Am 2 =7.69xl0~ 5 eV 2 and tan 2 0=0.45 [H). Dots represent 
the Borexino results from 7 Be and 8 B measurements. The 
error bars include also the theoretical uncertainty of the ex- 
pected flux from the Standard Solar Model BPS09(GS98). 



*^ = ^(^ e ).iV e .||(^). (7) 

where T e and E v are the electron and neutrino energies, 
and a x (x—e.fi-r) are the cross sections for elastic scatter- 
ing for different flavors. Tq—3 MeV is the energy threshold 
for scattered electrons, corresponding to a minimum neu- 
trino energy of E =3.2 MeV, N e is the number of target 
electrons, and d§ e /dE v is the differential 8 B solar neutrino 
flux [27j]. 

Using the above equation, we obtain P ee =0.29±0.10 at 
the mean energy of 8.9 MeV for 8 B neutrinos. Borexino is 
the first experiment to detect in real time, and in the same 
target, neutrinos in the low energy, vacuum dominated- and 
in the high energy, matter enhanced-regions. Borexino al- 
ready reported a survival probability of 0.56±0.10 for 7 Be 
neutrinos, at the energy of 0.862 MeV [9]. The distance 
between the two survival probabilities is 1.9 a. 

Future precision measurements of 7 Be and 8 B (and, possi- 
bly, pep) neutrinos in Borexino could provide an even more 
stringent test of the difference in P ee for low- and high- 
energy neutrinos predicted by the MSW-LMA theory: as- 
suming other 4 years of data taking, and to reduce the 
overall uncertainty on 7 Be neutrino rate at 5%, the distance 





Threshold 






[MeV] 


[10 6 cm" 2 s- 1 ] 


SuperKamiokaNDE 1 [3] 


5.0 


2.35±0.02±0.08 


SuperKamiokaNDE II [2] 


7.0 


2.38±0.05± 1 ^ 


SNO D 2 [4J 

SNO Salt Phase [25] 


5.0 
5.5 


9 -dq+0.24 +0.12 
^■jy_0.23 -0.12 

2.35±0.22±0.15 


SNO Prop. Counter [26J 


6.0 


i 77 +0.24+0.09 
L - 1 1 -0.21 -0.10 


Borexino 


3.0 


2.4±0.4±0.1 


Borexino 


5.0 


2.7±0.4±0.2 



Table VI: Results on B solar neutrino flux from elastic scat- 
tering, normalized under the assumption of the no-oscillation 
scenario reported by SuperKamiokaNDE, SNO, and Borex- 
ino. 

between the two survival probabilities can be improved at 3 
a. 
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